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ABSTRACT 


Under  this  program,  we  constructed  an  advanced  optical  phase  detection,  control  and 
analysis  system  for  space-time  interferometry  and  optical  communications.  The  system 
uses  scanning  interferometry  to  extract  the  phase  and  magnitude  of  serial  and  space-time 
fields.  Key  components  include  precision  mechanical  stages  and  controllers,  high  quality 
optical  elements,  serial  and  parallel  receiver  arrays,  broadband  sources,  and  advanced 
electronic  control,  data  acquisition  and  analysis  subsystems.  We  used  this  system  to 
extend  our  previous  work  on  space-time  image  capture  and  to  systematically  explore 
advanced  phase  characterization  techniques.  This  report  describes  our  results  in  using  this 
system  to  characterize  the  SNR  and  information  capacity  of  phase  sensitive 
interferometric  imaging  and  shows  experimental  evidence  of  super-resolved  scanning 
using  phase  sensitive  detection. 
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1.  INTRODUCTION 

Increasing  sophistication  in  the  control  and  detection  of  optical  fields  is  revolutionizing 
electro-optic  technology.  The  purpose  of  this  program  was  to  build  a  well  characterized 
interferometric  space-time  field  analysis  system  to  quantitatively  characterize  the 
information  capacity  of  emerging  detection  systems  for  applications  in  communications, 
microscopy,  and  data  storage.  A  space-time  system  uses  the  propagation  and  diffraction 
of  fields  to  transform  between  spatial  parameters  and  temporal  or  spectral  parameters. 
Interferometers  are  a  classic  example  of  systems  that  transform  spatial  parameters,  such 
as  mirror  position,  into  spectral  or  temporal  measurements  of  a  field.  Ultra-fast  pulsed 
lasers  and  optical  pulse  shapers  have  dramatically  expanded  the  variety  of  space-time 
systems  and  revitalized  the  need  for  advanced  serial  and  imaging,  linear  and  nonlinear 
interferometers  for  phase  sensitive  field  detection.  Over  the  past  several  years,  we  have 
developed  serial  and  imaging  cross-correlating  interferometers  for  phase  sensitive  field 
detection.  Instrumentation  provided  under  this  program  allowed  us  to  systematically 
characterize  the  capabilities  of  these  systems. 

The  space-time  analysis  system  built  under  this  program  has  components  in  four  key 
areas:  opto-mechnics,  data  acquisition,  broadband  sources  and  computational  analysis. 
Under  DURIP  support,  we  integrated  extra-ordinarily  precisely  controlled  and 
characterized  mechanical  positioning  systems  with  our  existing  interferometric  hardware 
and  we  dramatically  improved  our  data  acquisition  and  computational  analysis  hardware. 
These  improvements  allowed  US  to  quantitatively  evaluate  the  information  gathering 
potential  of  space-time  imaging  systems,  as  described  below.  We  also  acquired  a 
nonlinear  harmonic  generating  system  to  allow  us  to  shift  our  analysis  systems  to  key 
spectral  ranges  and  we  improved  our  existing  short  pulse  and  broad  band  source  systems. 

The  second  and  third  sections  of  this  report  detail  the  results  obtained  with  the  equipment 
supplied  under  this  grant.  Section  n  details  the  theory  of  interferometric  cross  correlation 
for  the  detection  of  space-time  fields  and  describes  experimental  demonstrations  of  noip 
and  information  limits  of  this  technique  probed  using  the  characterization  system  built 
under  this  program.  The  results  of  section  n  are  presented  in  the  refereed  literature  in 
“Noise  and  Information  in  Interferometric  Cross  Correlators,”  which  has  been  accepted 
for  publication  in  Applied  Optics.  The  authors  of  the  paper  are  Kent  B.  Hill  (now  at  the 
Naval  Research  Laboratory),  Scott  A.  Basinger  (now  at  the  Jet  Propulsion  Laboratory), 
and  Ronald  A,  Stack  and  David  J.  Brady  of  the  University  of  Illinois.  Section  m 
describes  the  results  of  particular  experimental  application  of  space-time  imaging  with 
interferometric  cross-correlators,  superresolved  optical  scaiming.  Section  HI  includes  a 
detailed  description  of  the  phase  and  amplitude  characterization  system  built  under  this 
program  and  describes  how  this  system  can  be  used  to  detect  subwavelength  shifts  in 
object  position  in  interference  microscopy.  The  results  of  section  HI  have  been  submitted 
for  publication  in  the  Journal  of  the  Optical  Society  of  America  A  as  “Superresolved 
optical  scanning  using  polychromatic  light,”  by  Scott  A.  Basinger,  Ronald  A.  Stack, 
Kent  B.  Hill,  and  David  J,  Brady. 
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II.  1  Noise  and  information  in  interferometric  cross  correlators 

Interferometric  cross  correlators  are  used  to  detect  the  amplitude  and  phase  of  opti¬ 
cal  signals  in  interferometric  microscopy ^  ultra-fast  signal  characterization^’^,  dispersive 
spectroscopy,  profilometry,  coherence  tomography,  and  time-gated  3-D  imaging.  The  noise 
properties  of  interferometric  spectrometers  have  previously  been  considered  in  the  context 
of  Fourier-transform  spectroscopy  (FTS),  where  the  focus  is  on  autocorrelation.  Developing 
applications,  such  as  interference  microscopy  and  coherence  tomography,  differ  from  FTS  in 
emphasizing  phase  sensitive  detection  and  cross  correlation.  The  most  important  difference 
between  auto  and  cross  correlation  measurements  is  that  while  the  time  bandwidth  product 
of  an  autocorrelation  is  essentially  one,  the  time  bandwidth  product  of  a  cross  correlation 
may  be  large.  The  additional  information  capacity  implied  by  this  difference  requires  that 
previous  studies  of  the  noise  characteristics  of  FTS  be  extended.  This  section  describes  phase 
sensitive  cross  correlation  and  evaluates  the  effect  of  various  noise  sources  both  theoretically 
and  experimentally. 

The  information  capacity  of  a  conventional  2-D  imaging  system  scales  linearly  with  the 
space-bandwidth  product,  meaning  that  as  one  increases  the  aperture  of  the  imaging  system 
one  expects  a  quadratic  increase  in  the  information  capacity.  In  principle,  the  information 
capacity  of  an  interferometric  imaging  system  should  be  linear  in  the  space-time  bandwidth 
product,  meaning  that  as  one  increases  the  temporal  scan  range  of  the  interferometer  one 
should  obtain  a  linear  increase  in  the  information  capacity.  Since  the  noise  in  the  cross 
correlation  increases  as  the  temporal  aperture  increases,  however,  the  information  detection 
capacity  of  a  given  correlator  ultimately  decreases  as  the  temporal  aperture  increases.  This 
meeins  that  there  exists  an  optimal  temporal  aperture  for  a  given  correlator  and  a  maximal 
information  capacity.  Analysis  of  these  quantities  is  a  key  result  of  this  section.  The  re¬ 
mainder  of  this  section  briefly  reviews  the  history  of  FTS  and  dispersive  FTS  and  motivates 
the  work  by  describing  some  applications.  Succeeding  sections  describe  the  process  of  in- 
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terferometric  cross  correlation,  analyze  noise  sources  and  information  capacity,  and  present 
experimental  results. 

Two  advantages  have  been  emphasized  in  comparing  FTS  with  other  approaches.  The 
Fellget  or  multiplex  advantage  is  that  the  signal-to-noise  ratio  for  a  measurement  of  M 
spectral  channels  may  be  improved  by  a  factor  of  y/M  relative  to  grating  spectrometry  due 
to  the  fact  that  all  spectral  elements  are  measured  simultaneously  in  FTS"^.  In  FTS  each 
spectral  element  is  measured  during  the  entire  measurement  time  T  rather  than  TfM  as 
is  done  in  a  grating  spectrometer.  The  Fellget  advantage  is  obtained  if  the  noise  detected 
during  the  measurement  does  not  increase  by  detecting  all  the  spectral  elements  simulta¬ 
neously.  The  Jacquinot  or  throughput  advantage  is  that  the  light-gathering  efficiency  is 
improved  compared  to  slit  spectrometers^.  A  thorough  analysis  of  FTS  was  presented  by 
Connes  in  1961^,  and  noise  levels  in  FTS  have  been  considered  by  several  workers  since 
then^^^^.  Cross  correlation  methods  analogous  to  the  techniques  discussed  in  this  section 
were  developed  later  in  the  1960s,  when  it  was  shown  that  a  sample  can  be  placed  in  one  arm 
of  an  interferometer  to  measure  its  wavelength  dependent  transmittance  or  reflectance^^’^^. 
Recent  reviews  of  this  technique,  known  as  dispersive  FTS,  have  been  given  by  Birch^"^  and 
Parker^^.  FTS  has  generally  been  used  in  the  far  infrared,  where  less  positioning  accuracy 
is  necessary  and  where  the  multiplex  advantage  is  significant.  In  the  near-IR  and  visible 
ranges,  FTS  has  no  fundamental  noise  advantage  over  grating  devices  but  still  maintains  the 
throughput  advantage.  The  throughput  advantage  is  of  particular  importance  in  imaging 
applications,  where  less  power  per  channel  is  available  and  where  high  spectral  sensitivity 
per  unit  analyzer  volume  is  vital.  As  translation  stage  technology  has  improved,  therefore, 
interferometric  techniques  for  visible  field  analysis  have  attracted  increasing  interest. 

Applications  of  interferometric  cross  correlators  include  profilometry  and  tomography. 
The  phase  of  the  field  scattered  from  an  uneven  or  rough  surface  provides  information  about 
the  depth  of  the  surface.  If  the  surface  depth  of  the  object  varies  slowly,  then  phase-shifting 
techniques  can  be  used  to  construct  a  map  of  the  surface  from  the  detected  signal^®’^^.  The 
resolution  of  these  measurements  is  often  limited  by  vibrations  that  cause  uncertainties  in 
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the  path  delay  between  the  reference  and  signal  beams.  Heterodyne  systems  can  be  used 
to  compensate  for  these  errors^®.  In  white-light  interferometry,  an  illuminating  source  with 
multiple  wavelengths  or  a  short  coherence  length  is  used  to  resolve  sharp  depth  variations 
on  a  surface^®”^®.  In  these  systems  the  location  of  the  central  fringe  in  a  cross  correlation 
determines  the  depth  of  the  object  from  which  light  is  reflected.  In  optical  coherence-domain 
reflectometry,  multiple  reflections  from  optical  components  are  detected^^  Similar  tech¬ 
niques  which  use  ultrashort  pulses  for  time-gated  imaging  of  objects  embedded  in  scattering 
media  have  been  developed^^’^®.  In  optical  coherence  tomography,  white-light  sources  can 
be  used  in  place  of  pulsed  sources,  and  the  short  coherence  length  acts  as  a  time  gate  for 
the  scattered  image^^®^ 

We  have  focused  on  two  emerging  applications  of  interferometric  cross  correlators:  repet¬ 
itive  ultra-fast  characterization  and  interferometric  microscopy.  Linear  interferometers  have 
not  generally  been  used  for  ultra-short  pulse  characterization  because  linear  autocorrelation 
yields  only  the  power  spectral  density,  rather  than  the  temporal  waveform  of  the  target  field. 
Nonlinear  techniques,  most  notably  frequency  resolved  optical  gating®^,  have  been  developed 
to  overcome  this  problem  and  determine  the  temporal  waveform  of  short  pulses.  Once  one 
has  determined  the  shape  of  a  short  pulse,  one  uses  the  pulse  to  probe  some  physical  situa¬ 
tion.  If  the  bandwidth  of  the  modified  probe  pulse  is  not  substantially  shifted,  the  temporal 
structure  of  the  modified  signal  can  be  determined  by  linear  cross  correlation  with  a  copy 
of  the  original  pulse.  This  technique  can  be  used  to  analyze  extraordinarily  complex  spaoe- 
time  fields,  as  demonstrated  recently  in  the  detection  of  holographically  encoded  space-time 
images^.  Interference  microscopy  expands  on  studies  of  white-light  interferometry  by  con¬ 
sidering  modulation  of  the  space-time  structure  of  fields  by  dispersive  scattering.  We  have 
shown  that  when  form  dispersion  is  considered,  a  3-D  interferometric  imaging  system  de¬ 
tects  more  information  about  scattering  objects  than  a  conventional  imaging  system^  In 
principle,  this  additional  information  should  lead  to  improvements  in  the  resolution  of  imag¬ 
ing  systems.  As  discussed  in  our  previous  work,  the  increase  in  the  detected  information  is 
highly  dependent  on  the  signal-to-noise  ratio  (SNR)  of  the  detection  system. 
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The  basic  process  of  detection  of  complex  optical  fields  by  interferometric  cross  correlators 
is  developed  in  Section  II. 2.  Signal-to-noise  analysis  for  practical  systems  and  a  discussion 
of  the  information  capacity  of  these  systems  are  presented  in  Section  II.3.  Experimental 
results  are  presented  which  demonstrate  the  detection  of  space-time  fields  and  their  noise 
characteristics  in  Section  II.4. 


II.2  INTERFEROMETRIC  CROSS  CORRELATION 

It  is  convenient  to  work  with  the  analytic  signal  representation  of  the  electromagnetic 
field.  The  analytic  signal  corresponding  to  a  real- valued  non-monochromatic  signal  (t) 
with  Fourier  spectrum  U^'''>{u)  is  defined  by  doubling  the  positive  frequency  components  and 
suppressing  the  negative  frequency  components^®’®'*.  Thus, 


u{t)  =  2  /  U^'^\u)exp{—j27:ut)du. 
Jo 


(1) 


The  Fourier  spectrum  of  the  analytic  signal  contains  no  negative  frequency  components. 

The  geometry  of  an  interferometric  cross  correlator  is  shown  in  Fig.  1.  A  reference 
field  whose  analytic  signal  is  given  by  u{x,y,t)  is  split  by  a  beamsplitter  in  a  Michelson 
interferometer.  In  one  arm  of  the  interferometer,  the  field  is  refiected  by  an  object  whose 
impulse  response  we  model  as  /io(ar,  y,  t).  The  signal  field  at  the  output  of  the  interferometer 
is  then  u{x,  y,  t)  *  ho  {x,  y,  t)  *hi{x,  y,  t) ,  where  ♦  denotes  convolution  and  hi  (x,  y,  t)  represents 
the  impulse  response  caused  by  diffractive  effects  of  the  rest  of  the  system.  In  the  second  arm 
of  the  interferometer,  the  reference  field  is  refiected  by  a  mirror  on  a  translation  stage  which 
controls  the  time  delay  through  this  path  and  to  the  detector.  If  both  paths  experience  the 
same  diffractive  effects  from  the  rest  of  the  system,  the  total  field  at  the  detector  is  given  by 


y'd{x,y,t)  =  v,h{x,y,t  -  r)  +  Uh{x,y,t)  *  ho{x,y,t),  (2) 

where  Uh  =  u*  hi  and  the  time  delay  r  reflects  the  difference  in  the  optical  path  length 
along  the  reference  and  signal  interferometer  arms.  The  detector  integrates  the  intensity  of 
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the  field  over  an  exposure  time  T,  and  the  output  signal  can  be  expressed  as 
fT  2 

I{x,y,T)=  /  \ud{x,y,t)\  dt 

Jo 

=  ^\uh{x,y,t)f  +  \ho{x,y,t)  *Uh{x,y,t)f 
+  ul{x,y,t  +  T)ho{x,y,t)  *Uh{x,y,t) 

+  Uh{x,y,t  +  T)  {ho{x,y,t)  *Uh{x,y,t)y]  dt.  (3) 


The  first  two  terms  on  the  right-hand  side  of  Eq.  (3)  are  constant  with  respect  to  r.  If  T 
is  long  enough,  then  the  last  two  terms  contain  the  temporal  cross  correlation  of  the  beams 
from  each  path  of  the  interferometer.  Since  interferometric  cross  correlation  is  a  first-order 
process,  the  contribution  of  the  incident  field  to  the  detected  temporal  signal  only  depends 
on  the  power  spectral  density  of  the  source.  If  the  input  field  is  spatially  uniform,  then  the 
last  two  terms  of  Eq.  (3)  can  be  written  as 


fT 

/  ul{t  +  T)ho{x,  y,t)  *  Uh{t)dt  -I-  ex. 

Jo 

=  f  dtu*f,{t  +  T)  f  dt'uh{t-t')ho{x,y,t')  +  c 

Jo  */— oo 

tT 

=  j  dt' ho{x,y,t')  dt  Uh{t  +  T)uh{t  —  t')  + 

=  [  dt' ho{x,y,t')Th{t' +  t)  +  C.C., 

J  —  OO 


.C. 


C.C. 


(4) 


where  c.c.  is  the  complex  conjugate  of  its  preceding  expression,  and  F/,  (r)  is  the  self-coherence 
function  of  the  optical  field  defined  as 


r(r)  —  f  u(t  -I-  T)u*{t)dt.  (5) 

It  can  be  seen  that  the  substitution  of  the  self-coherence  function  in  the  last  line  of  Eq.  (4) 
is  only  valid  if  T  is  much  longer  than  the  duration  of  the  impulse  response  and  the  width  of 
r(r).  According  to  the  Wiener-Khinchin  theorem,  the  self-coherence  function  can  also  be 
expressed  as  the  inverse  Fourier  transform  of  the  power  spectral  density  of  the  field: 

r(r)  =  f°°  \u{u)\'^ exp{-j2in4:)dj/.  (6) 

J—OO 
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Thus  a  temporally  incoherent  (white-light)  source  will  yield  the  same  intensity  as  a  coherent 
source  with  an  equivalent  power  spectral  density.  White-light  sources  can  have  much  wider 
bandwidths  and  are  much  less  expensive  than  mode-locked  laser  sources.  However,  it  is 
difficult  to  obtain  good  spatial  coherence  from  white-light  sources.  Spatial  coherence  is 
necessary  in  order  to  maintain  our  assumption  that  the  source  can  be  represented  as  a 
quasi-monochromatic  plane  wave. 

The  impulse  response  /io(x,  y,  t)  is  recovered  by  taking  the  Fourier  transform  of  /(x,  y,  r) 
with  respect  to  r.  For  u  >  0  and  assuming  that  Uh{t)  is  a  stationary  signal  and  that  T  is 
much  longer  than  the  duration  of  the  impulse  response  and  the  correlation  time  of  Uh  {t) ,  we 
obtain 

/  (x,  y,  u)  =  \Uh  Ho  (x,  y,  i/) ,  (^) 

where  we  make  use  of  the  fact  that  the  negative  frequency  components  of  the  analytic  signals 
are  zero.  Ho{x,y,u)  can  be  obtained  over  the  bandwidth  of  Uh{u)  by  dividing  Eq.  (7)  by 
l[/h(i/)p.  Given  Ho{x,  y,u),  we  can  reconstruct  a  bandlimited  version  of  the  object  impulse 
response  by  inverse  Fourier  transforming.  The  amount  of  information  that  can  be  extracted 
from  a  cross-correlation  measurement  depends  on  the  amount  of  noise  superimposed  on  the 
signal.  In  the  following  section,  we  consider  the  sources  of  noise  in  interferometric  cross 
correlators  and  determine  the  information  capacity  of  detected  signals. 


II.3  SNR  ANALYSIS  OF  INTERFEROMETRIC  DETECTORS 

In  this  section  we  consider  the  different  sources  of  noise  in  a  practical  interferometric  cross 
correlator  and  derive  SNR  limits  which  characterize  this  detection  technique.  We  begin  with 
a  qualitative  description  of  noise  sources  in  an  interferometric  cross  correlator.  Signal  noise 
is  caused  by  fluctuations  in  the  intensity  of  light  emitted  by  the  source  or  by  fluctuations  in 
the  transmissive  medium  between  the  source  and  the  detector.  In  an  interferometer,  signal 
noise  is  increased  by  mechanical  instabilities  which  cause  uncertainties  in  the  path  delay 
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between  the  interfering  beams.  Quantum  or  shot  noise  arises  in  the  photoelectric  detection 
of  light.  Thermal  noise  adds  to  the  signal  power  incident  on  the  detector  and  to  the  the 
signal  as  the  photocurrent  is  amplified  and  processed  by  analog  electronics.  Quantization 
error  is  caused  by  discrete  thresholding  of  the  signal  during  the  analog-to-digital  conversion 
process. 

The  relative  contribution  of  each  noise  source  depends  on  the  particular  system  and  its 
application.  Signal  noise  can  be  reduced  by  using  a  highly  stabilized  source  and  by  control¬ 
ling  the  environment  to  stabilize  fluctuations  in  the  transmissive  medium  in  the  interfer¬ 
ometer.  In  applications  where  the  visibility  of  the  interference  fringes  is  high,  the  intensity 
of  the  signal  varies  significantly  as  a  function  of  path  delay,  and  the  contribution  to  signal 
noise  from  inaccuracies  in  the  path  delay  increases.  This  effect  is  less  significant  in  the  IR 
region  where  the  intensity  of  the  fringes  varies  slowly  over  slight  changes  in  the  path  delay. 
However,  in  the  visible  region  of  the  spectrum  this  elBfect  must  be  considered.  Thermal  noise 
is  significant  when  the  photon  energies  of  the  optical  signal  are  comparable  to  or  less  than 
those  associated  with  thermal  radiation.  However,  detectors  of  visible  to  near-IR  radiation 
collect  very  little  thermal  radiation.  CCD  detectors  can  be  obtained  with  16-bit  analog-to- 
digital  converters,  and  the  thermal  noise  is  typically  less  than  one  bit  when  the  detectors  are 
cooled  to  low  temperatures  and  operated  at  narrow  electrical  bandwidths.  In  interferometric 
cross  correlators,  the  exposure  time  can  be  increased  to  decrease  the  bandwidth  and  reduce 
the  effects  of  thermal  noise.  With  a  dynamic  range  of  16  bits,  these  detectors  are  ultimately 
limited  by  shot  noise  when  detecting  visible  to  near-IR  radiation  at  moderate  to  high  signal 
intensities  and  narrow  bandwidths. 

We  now  proceed  with  a  quantitative  analysis  of  the  noise  in  an  interferometric  cross 
correlator.  To  recover  the  three-dimensional  signal  from  an  imaging  cross  correlator,  a  CCD 
array  is  used  to  detect  transverse  images  of  /(x,y,T)  at  evenly  spaced  intervals.  The  time 
delay  of  the  reference  beam  for  the  n*'*  sample  is  given  by  r„  =  nAr.  The  i,  m"*  pixel  in 
the  CCD  array  is  centered  at  {xi,y,n)  =  (I Ax,  m Ay)  and  has  area  A.  Assuming  that  the 
average  intensity  is  constant  over  the  area  of  each  pixel,  the  number  of  photons  collected  at 
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the  /,  m*'*  pixel  for  a  given  reference  delay  is  given  by 


Pl,m,n  —  l/m) '^n)» 

where  eo  is  the  free-space  permittivity  (in  c.g.s.  units),  c  is  the  speed  of  light,  wq  is  the 
center  frequency  of  the  field,  and  rj  is  the  quantum  efficiency  of  the  CCD.  For  notational 
convenience,  we  define 


7  = 


V 

2fia>o 


eocA. 


(9) 


In  practice,  shot  noise  contributes  to  the  number  of  detected  photons,  and  the  root  mean 
square  value  of  the  shot  noise  at  each  pixel  is  given  by  ■^pi,m,n'  Although  the  constant  terms 
in  Eq.  (3)  are  filtered  out  in  the  recovered  cross  correlation,  the  shot  noise  generated  by 
these  terms  is  evenly  distributed  over  the  entire  spectrum  of  the  detected  signal. 

The  noise  statistics  for  a  given  spatial  channel  in  our  system  are  independent  of  the 
spatial  location  of  the  channel.  Therefore,  the  x  and  y  dependencies  of  the  detected  signal 
are  neglected  in  the  following  analysis.  We  assume  that  for  a  given  path  delay  nAr,  the 
reference  field  at  the  plane  of  the  detector  can  be  written  as  u{t  -  nAr  +  Xn),  where  Xn 
is  the  positioning  error  in  the  reference  delay  for  the  sample.  The  number  of  photons 
detected  by  a  given  pixel  in  the  detected  image  can  then  be  written  as 


Pn  —  'yin  =  7  /  +  Xn)  +  ^(^)  * 

JQ 

where  a  is  a  scaling  constant  used  to  adjust  the  relative  intensities  of  the  signal  and  reference 
fields. 

The  signal  is  filtered  by  taking  discrete  Fourier  transforms  (DFT)  of  the  data  collected 
at  each  pixel  with  respect  to  r.  The  DFT  is  defined  by 

-Pn' =  SPn®’‘P(-^27rnn7Ar),  (H) 

n=:0 

where  Pn»  are  sampled  values  of  the  spectral  power  at  frequency  intervals  Ai/  =  1/iVAr. 
N  is  taken  to  be  even  for  simplicity.  In  calculating  the  statistics  of  we  consider  u{t) 
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as  a  stationary  random  process  and  Xn  as  a  random  variable.  We  assume  that  u{t)  is  a 
quasi-monochromatic  beam  and  that  the  variance  of  Xn  is  much  less  than  the  mean  optical 
period.  In  this  case,  u{t  —  nAr  +  Xn)  «  u  (t  —  nAr)  exp  (j27ri/oXn))  where  uq  =  a;o/27r.  We 
consider  first  the  expected  value  of  P„/.  Assuming  that  Xn  is  normally  distributed  with 
zero  mean  and  variance  and  that  the  expectation  values  of  the  reference  and  scattered 
intensities  are  independent  of  n,  we  find  for  n'  ^0  that  the  expected  value  of  Pn'  is 

/  N-\ 

(Pn')  =  (70:*  exp  (-27rVocr^)  exp  {j2nnn' I N) 

\  n=:0 

X  J  dt  j  du  exp{—j2nut)u“‘ {t  —  nAT)U  {i/)  Ho{u)^ 

/  N-l 

+  /7a  exp  (i27rnn7A’) 

\  n=0 

X  j  dt  j  du  exp  {—j27:ut)  u  {t  —  nAr)  U*  {—u)  H*  {—u)\  ,  (12) 


where  we  have  expressed  the  convolution  as  an  inverse  Fourier  transform.  Assuming  that  T 
is  large  enough  that  the  integration  range  over  t  can  be  considered  infinite,  we  can  perform 
this  integral  and  sum  to  obtain 


(Pn')  =  7a*  exp  (-2Tr^u^al)  J_^  du  exp  {N  -  1)  -  Atu^  ^ 

+  7a  exp  [—2'k^vIo^  j  du  exp  {N  —  1)  +  Atu 

sin(7r(Ari/  +  ^))  '  ' 


{Pn')  is  thus  the  convolution  of  the  target  spectrum  and  a  comb  function  and  is  periodic 
in  n'  with  period  N.  In  the  range  0  <  n'  <  N/2,  (P„/)  represents  (jU  {u)f'^Ho{u)  to  an 
approximate  resolution  of  1/NAt  as  long  as  the  period  of  the  comb  function  is  large  enough 


to  satisfy  the  Nyquist  criterion.  At  <  Atn  =  l/2i/raax,  where  Umax  is  the  largest  frequency 
at  which  ^[17  (t')l^)  Ho  {v)  is  nonzero. 
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Noise  arises  in  F„/  from  intensity  fluctuations,  shot  noise,  and  positioning  error  and 
instability.  We  estimate  the  impact  of  these  noise  sources  by  calculating  the  variance  of  P„/ . 
Since  P„'  is  a  complex  signal,  we  must  calculate  the  variances  of  the  real  and  imaginary 
parts  separately  and  their  joint  covariance.  It  can  be  shown,  however,  that  the  covariance 
of  the  real  and  imaginary  parts  vanishes  and  that  the  variances  of  the  real  and  imaginary 
parts  are  equal^'*.  Therefore,  it  is  sufficient  to  determine  the  quantity 

4  =  (l^’nf  -  |(/’n.)"l) 

=  (S  -  (S  {Pn.})"*)  +  (3  {P„.}*  -  (3  {P„'})'')  •  (14) 

Substitution  from  Eq.  (11)  yields 

4  =  E  {A  -  {?.>’>  =  E  4-  ('5) 

n=0  n=0 

The  total  variance  in  the  number  of  detected  photons  can  be  written  as^^ 

=  7(^n)  +7^1 

where  the  first  term  on  the  right-hand  side  is  contributed  by  shot  noise,  and  is  the  variance 
in  the  time- average  intensity  caused  by  intensity  fluctuations  and  positioning  errors. 

To  evaluate  Eq.  (15),  we  have  to  model  the  statistics  of  1^.  We  can  express  In  as 

/„  =  £?  +  exp  (-j27ri/oXn)  C'n  +  exp  {j2iruoXn)  C*,  (17) 

where 

E=\af  r\u{t)fdt+  r  \uit)*ho{t)fdt,  (18) 

Jo  Jo 

and 

Cn  =  a*  r  u*  {t  -  nAr)  u  (t)  *  ho  (t)  dt.  (19) 

Jo 

In  neglecting  the  significance  of  temporal  shifts  in  u{t),  we  have  assumed  above  that  E  is 
ergodic.  In  this  notation  the  variance  of  the  time  integrated  intensity  is 
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+  ((exp  (-i27ri/oXn)  Cn  +  exp  (j27ri^oXn)  C^f'^ 

-  (exp  {-j2'KV'oXn)  c„  +  exp  (j27rz^oXn) 

+  2{E  (exp  {-j2vuoXn)  Cn  +  exp  (j2xi^oXn)  Q)) 

-  2  (£;)  ((exp  i-j2iTUoXn)  Cn  +  exp  {j2iruQXn)  C*))  (20) 


The  first  term  in  Eq.  (20),  Cg,  is  independent  of  n.  For  quasi- monochromatic  illumi¬ 
nation,  Cn  and  are  weakly  modulated  harmonic  functions  of  n  with  no  dc  component. 
These  terms  thus  approximately  vanish  when  summed  over  n.  Substituting  Eq.  (20)  in 
Eqs.  (15)  and  (16),  thus  yields 

4  =  7^  N  (<7%  +  -(E))  -1-  ^  {(l^nl^)  -  l(exp  (-i27ri/oXn))|^  l(^n)|^}  .  (21) 

L  V  T'  /  n=0  -I 

Prom  Eq.  (13),  the  detected  signal  at  n'  =  NArUm  is  found  to  be 

l(Pn')l  =  (l^  (^tn)l^)  \HoM\ ,  (22) 

where  we  neglect  the  effects  of  aliasing.  If  necessary,  aliasing  could  be  reduced  digitally  by 
appropriate  apodization.  We  assume  that  Ho{y)  is  a  random  function  with  independent  real 
and  imaginary  parts.  The  mean  of  both  parts  is  assumed  to  be  zero  and  the  variances  are 
equal  to  .  The  SNR  for  individual  signal  components  (i.e.,  the  real  or  imaginary  part  of 
^fo(i/m)  for  a  particular  value  of  m)  is  the  ratio  of  |(P„'>|  /\/2  over  the  signal  statistics  to 


the  variance  of  the  real  part  of  P„<,  which  is  <tp/-\/2.  Thus, 

\a\  exp  (-27rVg4)  (\U  (i/m)|^) 


SNR  = 


At.  N  U  +  i(B))  +  E  {(|CnO  -  exp  (-2x^^<7=)  |(C„>|^} 


For  an  unpolarized  thermal  source,  the  variance  in  the  signal  intensity  is  given  by^ 


4  «  (24) 

where  Tc  is  the  coherence  time  of  the  signal.  For  a  white-light  source  T  >  Tc  and  this  term 
is  negligible  in  comparison  with  the  shot  noise  {E)/'y.  We  approximate  the  last  noise  term 
by  assuming  (|C„p)  «  |(C'n)P  and  writing 
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(25) 


Parseval’s  theorem  states  that  in  the  limit  NAt  -¥  oo, 


r  dt\c.\^=  r  du\T{Cn)\\  (26) 

7— oo  J—oo 

where  is  the  Fourier  transform  operator.  Using  the  convolution  theorem,  we  can  write 


:F{C„) 


(27) 


Therefore, 

/”  dt  =  lap  r  di,  \U{v)]*  |//.(!.)p .  (28) 

J-oo  J-oo 

Substituting  this  expression  into  Eq.  (25)  and  taking  the  average,  we  obtain 

E(|C„P>«J^(|t'(-^o)lf  4.  (29) 

where  ub  is  the  bandwidth  of  the  source. 

We  wish  to  calculate  the  SNR  as  a  function  of  the  various  signal  parameters  for  a  fixed 
total  detected  energy  per  exposure.  The  total  energy  that  can  be  detected  in  a  single  expo¬ 
sure  by  a  pixel  on  a  CCD  is  limited  by  its  well  capacity  P,,  or  the  number  of  photogenerated 
electrons  that  saturates  the  detected  signal.  If  the  duration  of  the  signal  is  long  compared 
to  Tc,  then  the  visibility  of  the  signal  fringes  is  weak  and  it  can  be  assumed  that  the  dc 
intensity  dominates  the  signal.  In  this  case, 

/E)  «  — .  (30) 

7 

Applying  Parseval’s  theorem  to  Eq.  (18)  we  obtain 

The  SNR  for  frequencies  within  the  bandwidth  of  the  source  can  be  estimated  by  setting 
i/„  =  i/fl.  Substituting  Eqs.  (29)-(31)  into  Eq.  (23),  we  find  that  the  SNR  is  maximized  when 
lap  =  and  we  obtain 
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SNR  = 


(32) 


_ exp  (-2TVg<7^) _ 

)■  ^  [l  _  exp  (-2xVgg»)l 

For  a  fixed  sampling  width  At  it  can  be  seen  that  increasing  N  causes  the  SNR  to 
decrease,  this  results  from  letting  more  noise  into  the  same  spectral  window.  However,  for  a 
fixed  measurement  duration  (the  product  of  AtN),  increasing  N  causes  the  expected  (x  \/jV 
increase  of  the  SNR.  One  may  also  improve  the  SNR  by  increasing  P,  or  the  dynamic  range 
of  the  CCD  camera  used.  To  increase  Pg ,  the  well  capacity  of  the  CCD  must  be  improved  or 
multiple  exposures  must  be  collected  for  each  measurement.  In  Fig.  2,  the  SNR  is  plotted  as 
a  function  of  the  variance  in  the  path  delay  to  demonstrate  the  effects  of  these  methods  for 
increasing  the  SNR.  To  calculate  the  SNR  plotted  in  Figs.  2-4,  the  source  was  assumed  to 
have  a  25  nm  bandwidth  centered  on  570  nm.  In  Fig.  2(a),  the  SNR  is  plotted  for  a  signal 
that  is  sampled  at  the  Nyquist  rate  with  AT  =  10^  and  P,  =  5  x  10®.  These  parameters  are 
typical  for  the  experimental  system  described  in  Section  II.4.  For  Fig.  2(b),  the  signal  is 
sampled  at  twice  the  Nyquist  rate  with  N  doubled  to  sample  the  same  signal  duration  and 
Pg  is  unchanged.  It  can  be  seen  that  the  SNR  increases  as  long  as  the  error  in  the  path  delay 
is  small.  This  corresponds  to  the  regime  in  which  shot  noise  dominates.  When  <7)^  is  large, 
noise  caused  by  the  positioning  error  dominates  and  the  SNR  does  not  increase  significantly. 
For  Fig.  2(c),  the  signal  is  sampled  at  the  Nyquist  rate  with  Pg  doubled  and  N  unchanged. 
In  this  case,  it  can  be  seen  that  the  SNR  improves  for  larger  values  of  In  both  methods 
for  improving  the  SNR,  the  total  energy  collected  during  the  measurement  was  doubled. 
This  accounts  for  a  decrease  in  the  shot  noise  since  shot  noise  scales  with  the  square  root 
of  the  number  of  photons  detected.  If  Pg  is  increased  by  taking  multiple  exposures,  it  is 
assumed  that  the  path  delay  remains  constant  during  all  of  the  exposures. 

In  Fig.  3,  the  SNR  is  plotted  as  a  function  of  N  for  two  values  of  Nyquist  sampling 
is  assumed  so  that  N  corresponds  to  the  signal  duration,  and  P,  =  5  x  10®.  When  the 
duration  of  the  signal  is  long,  shot  noise  dominates  due  to  the  large  dc  component  in  the 
detected  signal  and  the  positioning  error  has  little  effect  on  the  SNR.  The  positioning  error  is 
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more  significant  for  shorter  signals  because  the  amplitude  of  the  fringes  is  greater,  therefore, 
small  errors  in  the  path  delay  cause  significant  errors  in  the  detected  signal.  The  transition 
between  regimes  in  which  shot  noise  dominates  and  positioning  error  dominates  can  be  seen 
in  Fig.  4,  where  the  SNR  (normalized  to  its  maximum)  is  plotted  as  a  function  of  for 
several  values  of  N.  Nyquist  sampling  and  P,  =  5  x  10^  are  assumed.  It  can  be  seen  that 
for  signal  durations  longer  than  10^  periods  (approximately  100  ps),  the  positioning  error  is 
only  significant  when  it  is  greater  than  1%  of  the  wavelength.  For  shorter  signals,  smaller 
positioning  errors  become  much  more  significant. 

The  model  that  we  have  derived  for  the  SNR  in  interferometric  cross  correlators  provides 
an  understanding  of  the  limitations  of  these  systems.  Positioning  errors  and  mechanical 
instabilities  must  be  carefully  avoided  when  ultrafast  signals  are  detected.  To  minimize  the 
effects  of  shot  noise,  the  energy  collected  during  a  measurement  should  be  made  as  large  as 


possible. 

From  the  SNR  we  can  also  estimate  the  information  capacity  of  an  interferometric  cross 
correlator.  Let  SNRmi„  be  the  minimum  SNR  below  which  any  information  is  lost  in  the 
noise.  The  total  information  capacity  is  then 


rr  SNR 


(33) 


where  T  =  NAt  is  the  total  duration  of  the  measurement.  Figure  5  shows  the  information 
capacity  versus  T  for  Aq  =  566  nm,  <j^  =  6  nm,  Ps  =  500, 000,  ub  =  24  THz,  At  =  100  nm, 
and  SNRmin  =  80.  Note  that  the  information  capacity  has  a  maximum  at  a  certain  signal 
duration  and  decreases  as  the  length  of  the  signal  increases.  Thus,  increasing  the  information 
content  of  the  signal  may  actually  decrease  the  detected  information. 


II.4  EXPERIMENTAL  RESULTS 

To  demonstrate  the  detection  of  space-time  fields,  we  detected  the  field  reflected  from 
a  thick  scattering  sample  by  using  the  arrangement  of  Fig.  1.  A  beam  from  a  spatially 
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coherent  white-light  source  was  the  input  to  an  imaging  interferometric  cross  correlator  as 
described  in  Section  II. 2.  A  point  fluorescent  source  similar  to  that  developed  by  Liu  et  al.^® 
was  used  to  obtain  a  broad  bandwidth  and  good  spatial  coherence.  A  low  numerical  aperture 
lens  focused  the  image  scattered  from  the  object  with  high  depth  of  field  onto  a  scientific- 
grade,  cooled  CCD  with  16-bit  dynamic  range  (Princeton  Instruments)  at  the  output  of 
the  interferometer.  An  electro-motive  translator  (MPB)  with  external  feedback  from  an 
inductive  position  sensor  (Kaman  Instrumentation)  was  used  in  combination  with  a  linear 
translation  stage  (Aerotech)  to  control  the  path  delay  of  the  reference  beam.  Subwavelength 
steps  with  «  20-nm  absolute  resolution  were  obtained  over  40-^m  intervals,  and  the  linear 
stage  was  translated  by  40  nm  at  the  end  of  each  interval.  Neutral  density  filters  i^i,  i^2  und 
i^3  equalize  the  total  energy  collected  from  each  arm  of  the  interferometer  and  adjusted  the 
total  energy  incident  on  the  CCD.  By  placing  a  filter  in  each  arm  of  the  interferometer,  the 
dispersion  caused  by  the  filters  was  compensated.  Reflections  from  the  faces  of  the  filters  and 
the  sample  were  separated  from  the  image  by  slightly  tilting  these  surfaces.  A  PC  controlled 
the  experiment  and  stored  the  data  collected  from  the  CCD  during  the  cross  correlation. 

The  noise  characteristics  of  this  detection  system  were  determined  by  statistical  analysis 
of  a  series  of  autocorrelation  measurements.  A  total  of  99  autocorrelations  were  performed 
with  data  collected  from  121  pixels  in  each  measurement.  The  imaging  lens  was  removed 
during  these  measurements  to  yield  a  uniform  beam  intensity  over  the  area  of  these  pixels. 
The  SNR  ((P„*) /ap)  calculated  from  these  data  is  plotted  in  Fig.  6.  In  Fig.  6(a),  the  statistics 
were  calculated  using  the  data  from  each  of  the  121  pixels  in  a  single  measurement.  In  a  single 
measurement,  the  noise  due  to  positioning  errors  and  mechanical  instabilities  is  identical  for 
each  pixel.  Therefore,  only  the  noise  due  to  source  intensity  fluctuations,  quantum  noise, 
thermal  noise  and  electrical  noise  is  included  in  the  SNR  plotted  in  Fig.  6(a).  The  maximum 
SNR  in  this  plot  is  approximately  250.  In  Fig.  6(b),  the  noise  due  to  all  effects  was  included 
by  calculating  the  statistics  of  the  spectral  data  collected  from  a  single  pixel  over  all  of  the 
measurements.  In  this  case,  the  SNR  was  reduced  to  about  200. 

It  appears  that  the  primary  source  of  noise  in  our  experiments  was  shot  noise.  The 
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average  number  of  photoelectrons  detected  in  each  data  point  was  (pn)  ^  2.2  x  10^,  and  the 
number  of  samples  taken  was  N  =  309.  By  substituting  our  experimental  parameters  into 
Eq.  (32),  the  theoretical  value  of  the  SNR  is  found  to  be  approximately  900  for  no  positioning 
error.  With  a  positioning  error  of  1%  of  the  wavelength,  the  theoretical  SNR  drops  to  about 
180.  To  improve  the  SNR,  it  is  necessary  to  increase  (p„);  however,  the  saturation  level 
for  our  CCD  is  5  x  10®  photoelectrons.  Averaging  several  frames  collected  by  the  CCD 
would  reduce  the  effects  of  shot  noise,  but  the  time  required  for  experiments  would  increase 
proportionately.  In  this  case,  positioning  errors  would  have  a  greater  effect  unless  the  path 
delay  can  be  held  constant  during  the  time  required  to  average  several  frames. 

Systematic  errors  in  the  data  collected  during  these  experiments  were  eliminated  by 
manipulating  the  data  to  compensate  for  the  errors.  Each  measurement  lasted  approximately 
5  min,  and  the  total  time  required  to  perform  99  measurements  was  almost  8  h.  Because  the 
beam  intensity  was  not  perfectly  uniform  over  the  area  of  the  pixels  used  to  collect  the  data, 
and  because  the  average  intensity  of  the  source  varied  from  one  measurement  to  the  next, 
the  data  from  each  pixel  in  each  measurement  was  normalized  to  have  zero  mean  and  unity 
standard  deviation.  Furthermore,  because  slowly  varying  temperature  fluctuations  during 
the  experiments  caused  the  optical  path  difference  to  change  from  one  measurement  to  the 
next,  the  spectrum  in  each  measurement  was  multiplied  by  a  phase  factor  to  compensate 
for  the  corresponding  time  shift  in  the  autocorrelation. 

To  experimentally  verify  the  theoretical  analysis  of  the  SNR  as  a  function  of  positioning 
error  would  require  several  thousand  measurements,  with  intentional  positioning  noise  added 
to  the  reference  mirror  sampling  positions.  However,  this  may  not  yield  satisfactory  results 
due  to  the  inherent  positioning  error  already  contained  in  the  system.  As  an  alternative,  we 
took  one  large  set  of  highly  oversampled  data  by  moving  the  reference  mirror  in  step  sizes 
of  20  nm.  Again,  we  collected  data  from  121  pixels.  In  addition,  we  added  a  scheme  to 
interferometrically  detect  the  optical  path  difference  between  the  reference  mirror  and  the 
object  mirror  to  less  than  1  nm  over  the  full  range  of  motion  of  the  reference  mirror  (4  inches) . 
This  scheme  uses  a  reference  laser  beam  (Spectra  Physics  frequency  stabilized  HeNe)  which 
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is  placed  slightly  above  and  parallel  to  the  signal  beam.  The  spatial  fringe  pattern  produced 
from  this  HeNe  is  then  captured  with  the  CCD  array  and  Fourier  transformed  to  extract  the 
phase.  By  maintaining  a  careful  history  of  the  HeNe  phase,  the  absolute  reference  mirror 
position  can  be  calculated.  From  this  real  oversampled  data  set,  we  then  simulated  data  sets 
with  various  positioning  errors.  Gaussian  random  variables,  with  a  standard  deviation  of 
the  desired  positioning  error,  were  added  to  uniformly  sampled  positions  and  the  intensity 
was  then  interpolated  at  these  positions  from  the  real  data  set.  Figure  7  shows  the  results 
for  three  different  values  of  the  total  number  of  samples  taken  N.  In  each  case,  the  total 
length  of  each  measurement  T  was  held  constant,  so  that  a  greater  value  of  N  corresponds 
to  a  higher  sampling  density,  which  is  why  the  SNR  increases  with  N.  As  a,  comparison,  the 
dashed  lines  show  the  same  corresponding  calculation  for  the  SNR  versus  positioning  error 
as  determined  by  Eq.  32. 

For  the  measurement  of  3-D  cross  correlations,  a  scattering  sample  was  made  by  adding 
5-/xm  alumina  polishing  powder  to  a  transparent  epoxy  (Epo-Tek  301).  The  epoxy  was  then 
poured  between  two  microscope  cover  slides  and  allowed  to  harden.  The  thickness  of  the 
sample  was  4  mm,  and  the  cover  slides  were  150-/xm  thick.  A  cross  correlation  of  the  field 
scattered  from  this  sample  was  measured  by  collecting  242  /xm  x  242  /xm  images  on  the  CCD 
(11  X  11  pixels)  at  a  rate  of  approximately  four  samples  per  fringe  as  the  reference  delay 
T  was  increased.  Because  of  the  difficulty  in  plotting  3-D  data,  density  plots  of  only  three 
2-D  slices  of  the  data  are  shown  in  Fig.  8.  The  drawing  at  the  top  of  the  figure  shows  the 
locations  of  the  2-D  slices  that  are  plotted  below  it.  For  clarity,  and  to  reduce  the  amount 
of  data  to  display,  only  the  first  5.6  ps  of  the  scattered  field  was  detected.  A  total  of  10,752 
images  were  collected  over  this  range.  The  DFT  of  the  data  collected  from  each  pixel  was 
calculated,  and  the  inverse  DFT  of  350  points  centered  at  the  fringe  frequency  provided  the 
filtered  data  shown.  At  r  «  0.2  ps,  the  light  scattered  from  the  first  surface  of  the  sample 
appears,  and  the  tilt  of  this  surface  can  be  seen.  No  further  scattering  can  be  seen  until 
r  >  1.52  ps,  which  corresponds  to  the  thickness  of  the  glass  cover  slide.  The  data  collected 
from  one  pixel  in  the  image  is  shown  in  Fig.  9(a).  The  DFT  and  the  amplitude  of  the  filtered 
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data  for  this  pixel  are  plotted  in  Figs.  9(a)  and  (b),  respectively. 


II.5  CONCLUSION 

We  have  shown  that  interferometric  cross  correlators  can  be  used  to  detect  complex 
space-time  fields.  Because  an  interferometric  cross  correlation  depends  only  on  the  power 
spectral  density  of  the  source,  the  impulse  response  of  an  object  can  be  determined  using 
simple  white-light  sources.  To  recover  the  phase  of  the  signal,  sampling  of  the  cross  correla¬ 
tion  at  the  Nyquist  rate  must  be  performed.  In  the  visible  region,  the  main  sources  of  noise 
in  interferometric  cross  correlators  are  shot  noise  and  inaccuracies  in  the  time  delay  of  the 
reference  field.  The  SNR  in  these  systems  decreases  exponentially  with  the  variance  of  the 
error  in  the  path  delay  and  is  scaled  by  a  factor  of  1/  ^/N ,  where  N  is  proportional  to  the 
length  of  the  signal.  We  have  developed  hardware  for  an  interferometric  cross-correlation 
system,  and  the  SNR  achieved  in  this  system  is  approximately  200. 
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FIGURES 


Fig.  1.  Experimental  setup  for  detecting  the  impulse  response  of  a  thick  scattering  object. 

Fig.  2.  SNR  vs.  path  delay  error,  (a)  N  =  10“*,  At  =  Arpf,  and  =  5  x  10®.  (b)  iV  =  2  x  10“*, 

At  =  ATff/2,  and  P,  =  5  X  10®.  (c)  =  10^  At  =  At/v,  and  P,  =  10®. 

Fig.  3.  SNR  vs.  signal  duration  for  two  values  of  positioning  error. 

Fig.  4.  Normalized  SNR  vs.  positioning  error  for  different  signal  durations. 

Fig.  5.  Information  capacity  as  a  function  of  signal  duration  for  Aq  =  566  nm,  <7^^  =  6  nm, 
Pg  =  500, 000,  VB  =  24  THz,  At  =  100  nm,  aind  SNRmin  =  80. 

Fig.  6.  SNR  for  experimental  cross-correlation  system,  (a)  SNR  without  positioning  and  me¬ 
chanical  errors,  (b)  Total  SNR.  The  dashed  line  shows  the  shape  of  the  mean  spectrum  with 
arbitrary  scale. 

Fig.  7.  Experimentally  determined  SNR  versus  positioning  error  for  three  different  sampling 
densities  and  the  corresponding  theoretical  calculations  (dashed  lines). 

Fig.  8.  Density  plots  of  the  space— time  field  scattered  from  a  thick  scattering  sample. 

Fig.  9.  Cross-correlation  data  from  one  pixel  in  the  space— time  field  scattered  from  a  thick 
scattering  sample,  (a)  Raw  data,  (b)  DFT  of  raw  data,  (c)  Amplitude  of  the  filtered  cross 
correlation. 
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.  1.  Experimental  setup  for  detecting  the  impulse  response  of  a  thick  scattering  object. 


Fig.  4.  Normalized  SNR  vs.  positioning  error  for  different  signal  durations. 


Fig.  5.  Information  capacity  as  a  function  of  signal  duration  for  Aq  —  566  nm,  <7-^ 
Ps  =  500, 000,  i/B  =  24  THz,  At  =  100  mn,  and  SNRmin  =  80. 
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Fig.  6.  SNR  for  experimental  cross-correlation  system,  (a)  SNR  without  positioning  and  me¬ 
chanical  errors,  (b)  Total  SNR.  The  dashed  line  shows  the  shape  of  the  mean  spectrum  with 
arbitrary  scale. 
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Fig.  7.  Experimentally  determined  SNR  versus  positioning  error  for  three  different  sampling 
densities  and  the  corresponding  theoretical  calculations  (dashed  lines). 
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Fig.  8.  Density  plots  of  the  space-time  field  scattered  from  a  thick  scattering  sample. 
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The  concept  of  “resolution”  is  based  on  the  observation  that  a  given  system  cannot 
distinguish  between  different  objects  if  the  differences  are  on  too  small  of  a  spatial  scale. 
Translation  is  among  the  sin^)lest  techniques  for  creating  different  objects.  A  classical 
imaging  system  cannot  distinguish  between  an  object  and  the  same  object  shifted  by  an 
amniint  less  than  the  classical  resolution  limit  of  one-half  the  center  wavelength  of  the 
probing  light.  The  resolution  limit  of  an  interferometric  imaging  system  is  discussed 
below.  The  goal  of  this  section  is  to  show  that  the  resolution  limit  of  an  interferometric 
imaging  system  using  polychromatic  light  can  exceed  that  of  a  classical  imaging  system. 
We  do  tto  by  demonstrating  that  our  interferometric  imaging  system  can  detect  much 
smaller  transverse  shifts  in  object  position  than  its  classical  counterpart  could. 

Our  expectation  that  an  interferometric  imaging  system  should  yield  higher 
resolution  can  be  understood  by  considering  how  information  is  encoded  in  imaging  and 
tomographic  systems.  The  “object”  in  such  systems  is  a  passive  spatial  structure  which 
scatters  probe  or  ambient  light  into  the  imaging  system.  The  anqilitude  and  phase  of  the 
scattered  light  vary  as  functions  of  spatial  structure  and  con^osition  of  the  object.  The 
system  collects  this  scattered  light  and  uses  it  to  reconstruct  an  “image,”  which  is  a  map  of 
the  spatial  and  spectral  distribution  of  the  scattered  field  in  the  neighborhood  of  the  object. 
The  structure  of  the  object  is  inferred  by  analogy  with  the  image.  In  a  classical  imaging 
system,  the  relative  phase  of  different  tenqioral  frequency  components  is  lost  on  detection. 
The  relative  phase  is  particularly  in^ortant  in  determining  the  effects  of  form  dispersion 
which  can  be  used  for  subwavelength  detection.  A  high  spatial  frequency  structure  in  an 
object,  while  not  directly  observable  in  its  effect  on  the  intensity  of  the  scattered  field,  can 
shift  the  phase  of  spectral  and  spatial  con^onents  in  the  scattered  field.  These  phase  shifts 
can  be  detected  by  considering  the  relative  phases  between  the  different  spectral 
con^onents  in  the  scattered  field.  Interferometric  detection  is  necessary  to  detect  these 
relative  phases. 

Since  the  scattering  processes  associated  with  imaging  are  linear,  it  is  not 
necessary  to  detect  the  temporal  characteristics  of  the  scattered  field  to  obtain  the  relative 
phases  between  different  frequency  con^nents.  These  relative  phases  can  be  determined 
by  measuring  the  impulse  response  of  the  scattering  process.  By  “impulse  response”  we 
mean  the  field  scattered  when  the  object  is  illuminated  by  an  impulse  in  the  spatial 
frequency  and  temporal  fi^uency  domains,  e.g.  a  plane  wave  temporal  inqiulse.  In 
practice,  it  is  not  possible  to  generate  or  detect  the  complete  in^ulse  response,  but  a 
ff.Tnpnrally  and  spatially  band-hmited  version  can  be  detected  interferometrically  using  a 
broadband  illumination  source.  The  use  of  a  Michelson  interferometer  to  determine  the 
impulse  response  of  a  scattering  object  is  discussed  in  Ref.  [5].  Using  this  technique,  one 
should  be  able  to  increase  the  amount  of  information  detected  from  an  object  by  increasing 
the  temporal  extent  of  the  impulse  response  measured.  This,  in  turn,  increases  the 
resolution  of  the  system.  However,  there  is  a  limit  to  the  amount  of  information  that  one 
can  detect  interferometrically.  As  discussed  in  Ref.  [1],  the  signal-to-noise  ratio  (SNR)  of 
an  interferometric  system  limits  the  amount  of  information  one  can  extract  from  a 
measured  impulse  response.  Thus,  the  resolution  limit  of  an  interferometric  imaging 
system  is  bounded  by  the  SNR  of  the  system.  The  SNR  of  the  system  described  here  is 
characterized  in  Ref.  [5] 
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In  this  section,  we  describe  experimental  results  confirming  that  interferometric 
imaging  of  a  scattered  polychromatic  optical  field  can  be  used  to  detect  superresolved 
spatial  structures.  Section  2  discusses  the  data  analysis  protocol  underlying  our 
experiments.  As  mentioned  above,  inversion  of  3-D  data  collected  by  our  interferometric 
imaging  system  to  obtain  images  is  beyond  our  current  capability.  The  data  analysis 
protocol  is  the  model  we  use  to  distinguish  scattering  data  collected  for  different  object 
positions.  Section  3  describes  the  experimental  system  and  general  procedure  used  to 
collect  scattering  data.  Finally,  Section  4  presents  the  experimental  results. 


m.2.  ANALYSIS  PROTOCOL 

In  this  section,  the  process  of  analyzing  data  collected  with  our  interferometric  cross- 
correlator  is  discussed.  The  purpose  of  our  analysis  is  to  discriminate  between  data  sets 
collected  at  several  different  object  positions.  To  gain  an  understanding  of  the  principles 
and  motivation  behind  our  analysis,  several  aspects  of  our  system  and  data  collection 
procedure  must  be  explained.  First,  the  motivation  as  to  the  type  of  data  taken  is 
discussed  and  then  the  general  method  of  data  acquisition  is  described.  Finally,  we  show 
how  the  data  is  processed  to  extract  features  and  measure  object  position. 

As  stated  in  Section  1,  “classical”  imaging  systems  measure  the  time-averaged 
intensity  of  an  optical  field,  thus  losing  degrees  of  fireedom  associated  with  the  ten5)oral 
structure  of  the  field.  Ref.  [1]  argues  that  the  temporal  degrees  of  fireedom  that  are 
encoded  in  the  phase  spectrum  of  the  scattered  optical  field  need  to  be  recovered  in  order 
to  resolve  a  subwavelength  spatial  structure.  While  the  information  received  by 
independent  spatial  and/or  spectral  channels  is  not  conq)letely  independent,  some 
information  associated  with  a  subwavelength  spatial  structure  is  recovered  by  analyzing 
the  enseihble  of  data  captured  from  the  phase  spectrum  of  the  scattered  optical  field  over  a 
spatially  separated  array  of  detector  channels.  There  is  no  direct  linear  mapping  between 
this  information  and  the  spatial  structure  of  an  object’s  surface.  Since  the  information 
received  in  different  channels  is  not  independent,  doubling  the  nuniber  of  spatial  channels 
or  doubling  the  number  of  spectral  conqjonents  does  not  necessarily  double  the 
information  recovered.  The  inter-dependence  of  the  data  and  the  lack  of  a  linear  mapping 
malcft  image  construction  fi’om  polychromatic  data  extremely  challenging. 

For  our  analysis,  we  simplify  the  problem  considerably  by  attempting  not  to  form 
an  image  of  an  object,  but  only  to  distinguish  between  several  different  object  positions. 
This  is  a  sin^ler  task,  since  we  can  employ  statistical  methods  on  this  data  to  find 
correlations  among  different  object  positions.  For  example,  assume  that  one  has  several 
sets  of  phase  spectrum  data  from  one  position,  and  several  more  sets  from  another 
position.  If  we  can  select  one  set  from  the  ensemble  and  determine  which  object  position 
it  came  fix)m,  we  have  succeeded  in  demonstrating  that  our  system  is  capable  of 
subwavelength  resolution. 

The  statistif  fll  correlations  we  seek  rely  on  data  sets  containing  the  relative  phases 
between  different  spectral  conponents  of  the  scattered  field  over  many  spatial  channels. 
With  an  interferometer  we  perform  a  cross-correlation  that  results  in  an  intensity  at  the 
detector  for  the  nu,  sanple  with  the  following  form. 
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(1) 


J'T 

(x,  y,t-n8T)  +  u^ix,y,  t)*  h„  (x,  y,t)dt. 

The  function  Uhix,  y,  t)  is  the  in^ulse  response  of  the  source,  and  ho(x,  y,  t)  is  the  impulse 
response  of  the  object.  The  variable  St  is  the  temporal  sanpling  width  and  **’  represents 
convolution.  Note  that  x  and  y  are  also  discrete  in  the  sense  that  they  correspond  to 
individual  spatial  channels  (pbtels)  on  the  detector.  This  set  of  data  hix,  y),  n  =  l,...,N , 
where  N  is  the  total  number  of  san^les,  is  not  particularly  useful.  However,  in  Ref.  [5]  it  is 
shown  that  a  Fourier  transformation  of  this  data  set  can  be  used  to  obtain  the  transfer 
function  Ho(x,  y,  v)  of  the  object.  More  specifically,  performing  a  fast  Fourier  transform 
(FFT)  on  I„(x,  y)  with  respect  to  n,  results  in  the  object  transfer  function  times  the  power 
spectrum  of  the  source. 


In(x,y)  = 


(2) 


Given  that  the  power  spectrum  of  the  source  is  a  known  quantity,  it  can  be  divided  out 
leaving  only  the  transfer  function.  Since  the  object  transfer  function  is  the  same  as  the 
phase  and  anplitude  spectrum  of  the  object,  a  Fourier  transformation  of  the  cross¬ 
correlation  data  will  give  the  relative  phase  information  form  dispersion  necessitates.  Note 
that  the  object  transfer  function  times  the  power  spectrum  of  the  source  will  be  referred  to 
fi’om  this  point  as  the  transfer  function  since  the  power  spectrum  term  does  not  enter  into 
our  statistical  analysis. 

An  imaging  interferometric  cross-correlator  consists  of  the  following:  a  point 
white-light  source  which  is  collimated  and  used  to  probe  or  image  the  object  in  question,  a 
Michelson  interferometer  which  allows  us  to  interfere  a  reference  drawn  fi'om  the  probe 
light  with  light  reflected  fi’om  the  object,  and  a  detector  (charged  coupled  device  or  CCD) 
array  which  captures  the  interference  pattern.  An  interferometer  performs  the  cross¬ 
correlation  by  adjusting  the  distance  traveled  by  the  reference  beam.  The  difference 
between  distance  traveled  by  the  reference  beam  and  the  distance  traveled  by  the  object 
beam  is  called  the  optical  path  difference  (OPD).  Scanning  the  OPD  throughout  the  entire 
range  where  an  interference  pattern  can  be  detected  while  sanq)ling  the  intensity  over 
several  spatial  channels  results  in  the  set  of  cross-correlation  data.  The  range  of  the  OPD 
is  chosen  to  cover  the  temporal  extent  of  the  band-limited  object  signal,  while  the  spacing 
is  chosen  such  that  the  san^)ling  rate  is  close  to  the  Nyquist  limit  of  the  temporal 
interference  fiinges.  This  cross-correlation  data  is  collected  numerous  times  for  several 
different  object  positions. 

Our  statistical  analysis  consists  of  reducing  the  large  amount  of  data  collected  to  a 
small  set  of  averaged  parameters,  a  set  for  each  object  position.  Comparing  the  absolute 
value  of  the  difference  between  the  average  set  of  parameters  for  a  particular  position  j  to 
an  individual  set  of  parameters  for  a  particular  position  i  yields  a  number  which  reveals 
how  well  the  data  at  position  i  matches  the  ensemble  of  data  sets  taken  at  a  position  j. 
Determining  when  this  number  is  a  minimum  will  indicates  at  which  position  this  data  was 
taken  (i.e.  when  i  =  j.)  A  detailed  description  of  this  analysis  follows. 

The  cross-correlation  data  is  interpolated  to  evenly  sampled  points  using 
accurately  measured  OPD  values.  This  is  done  since  our  system’s  ability  to  accurately 
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measure  the  OPD  is  better  than  its  ability  to  control  the  OPD  or  to  position  the  reference 
mirror.  Interpolating  to  evenly  spaced  sample  points  greatly  simplifies  the  task  of 
performing  a  Fourier  transform  on  the  data,  since  its  allows  an  FFT  to  be  used.  Carrying 
out  an  FFT  results  in  a  transfer  function  for  the  object  which  is  a  complex-valued  set  of 
data.  The  data  is  square-window  filtered  to  keep  only  those  fi^equency  components  which 
are  present  in  the  source  spectrum.  This  reduces  both  noise  and  the  amount  of  data  to 
process  in  following  steps.  The  data  is  now  arranged  into  an  m  vector  f,  where  m  is  the 
total  number  of  frequency  components  kept  times  the  number  of  spatial  channels  captured 
per  cross-correlation.  These  data  sets,  f,  ,i  =  1,. . .  ,n ,  where  n  is  the  total  number  of  cross¬ 
correlations  captured,  are  arranged  into  an  m  x  n  matrix  F: 

. f,].  (3) 

The  columns  of  F  now  contain  the  transfer  functions  of  the  object.  To  use  statistical 
correlations,  some  of  the  columns  contain  cross-correlations  collected  with  an  object  at  a 
certain  position,  call  it  position  1,  other  columns  contain  cross-correlations  collected  at 
position  2,  and  so  on. 

To  make  correlations  between  the  columns  F(f, )  easier  to  find,  we  need  to  reduce 
the  amount  of  data  to  analyze  without  losing  any  information.  One  useful  method  is 
principle  component  analysis.  First,  we  perform  a  singular  value  decon5)osition  on  F  to 
find  an  orthonormal  set  of  basis  vectors  U/,i  =  for  each  f,.  These  basis  vectors 

are  arranged  in  descending  order  according  to  their  corresponding  singular  values. 
Principal  conq)onent  analysis  tells  us  that  only  the  first  few  basis  vectors  U,  are  necessary 
to  represent  each  f,  to  within  a  certain  tolerance.  Therefore,  we  choose  to  truncate  this 
expansion  to  a  length  of  /i^  sufficient  for  the  accuracy  of  this  analysis.  Now  the  original 
vectors  f  j  can  be  expressed  in  terms  of  the  basis  vectors  U,  as 


f.=S“uU/> 


(4) 


where  a,  j  =  f/U*  are  the  expansion  coefficients. 

Calculating  the  expansion  coefficients  for  each  f,  gives  us  the  reduced 

set  of  data  desired.  Each  cross-correlation’s  transfer  function  is  now  represented  by  only 
five  or  six  expansion  coefficients.  Since  each  position  of  the  object  has  multiple  sets  of 
cross-correlations,  the  expansion  coefficients  for  identical  object  positions  are  averaged 
together.  We  define  the  average  of  the  expansion  coefficients  for  each  position  k  as 

•  One  way  to  quantitatively  find  statistical  correlations  among  the  recorded 
cross-correlations  is  to  con^are  the  expansion  coefficients  for  cross¬ 
correlation  with  the  average  expansion  coefficients  for  each  object  position: 
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(5) 


I 

Xi.*  ~  j  * 
y=i 

The  object  position  k  at  which  an  individual  cross-correlation  i  was  taken  should  provide  a 
minimum  for  that  particular  value  of  Yi,k  •  P®*"  example,  if  the  14th  cross-correlation  f  ,4  is 

from  object  position  3,  then  yj4  3  should  be  a  minimum  since  the  expansion  coefficients 
should  be  similar  to  the  average  expansion  coefficients  for  that  object  position.  For 
{Tm*}  ’  expect  larger  values.  In  this  way,  once  we  have  collected  several 

cross-correlation  from  different  object  positions,  we  can  identify  which  object  position 
subsequent  cross-correlation  are  collected  from  by  looking  for  a  minimum  Yi,k  • 


m.3.  EXPERIMENTAL  SETUP  AND  PROCEDURE 

We  have  experimentally  verified  the  validity  of  the  analysis  presented  above  and  in  Ref.  [1] 
by  building  the  interferometric  imaging  system  outlined  in  Fig.  1.  As  explained  in 
Section  2,  an  interferometer  is  used  for  performing  cross-correlations  between  a  reference 
and  object  signal  In  turn,  this  cross-correlation  is  used  to  find  the  transfer  function  for  the 
object.  Here,  we  discuss  the  experimental  setup  of  our  Michelson  interferometer  which 
we  use  to  find  the  transfer  function  produce  by  an  object  illuminated  by  a  spatially 
coherent  polychromatic  or  broadband  source. 

A  point  fluorescent  source  similar  to  that  developed  by  Liu  et  al.’^  is  used  to  obtain 
a  broad  spectral  bandwidth  with  acceptable  spatial  coherence.  The  source  is  a  dye-jet 
illuminated  by  the  multi-spectral,  2.5  W  output  of  an  Argon-ion  laser.  A  sapphire  nozzle 
in  the  dye-jet  provides  a  laminar  dye  stream  400  jim  thick.  The  dye  used  is 
Rhodamine  590  Chloride  (Exciton),  also  known  as  Rhodamine  6G.  The  dye  is  dissolved 
in  ethylene  glycol.  The  laser  is  focused  in  the  center  of  the  dye-stream  by  a  25  x 
microscope  objective  to  a  spot  size  of  1.6  pm.  The  idea  is  to  focus  the  laser  beam  to  the 
smallest  spot  possible  to  simulate  the  infinite  spatial  coherence  from  an  ideal  point  source. 
T  ight  from  the  fluorescing  dye  is  re-collimated  by  the  same  microscope  objective  and  then 
filtered  by  a  band-pass  absorption  filter  (Oriel  51302)  to  filter  out  the  original  laser  source. 
The  output  is  a  spatially  coherent,  broad-band  beam  having  a  spectrum  with  a  center 
wavelength  of  566  nm  and  a  FWHM  of  approximately  30  nm  The  difference  between  this 
light  source  and  the  one  modeled  in  our  theoretical  analysis  (Ref.  [1])  is  that  this  source 
has  a  continuous  spectrum  instead  of  several  distinct,  sharp  frequencies.  This  results  in  a 
lower  overaD  SNR  for  our  experimental  system,  but  provides  us  with  a  range  of 
fiequencies  with  which  to  probe  our  object. 

The  Michelson  interferometer  splits  the  collimated  beam  produced  by  the  dye 
source  into  an  object  beam  and  reference  beam  This  is  done  using  a  1-inch  beamsplitter 
cube  (OFR  SCL-25-VIR)  which  is  anti-reflection  (AR)  coated  to  pass  light  in  the 
550-780  nm  range.  Beyond  the  beamsplitter,  neutral  density  filters  are  placed  in  each  path 
to  equalize  the  total  energy  reflected  from  the  reference  mirror  and  the  object.  The  reason 
for  placing  a  filter  in  each  arm  of  the  interferometer  is  to  neutralize  the  dispersion  an 
individual  filter  in  just  one  arm  would  cause.  Reflections  from  the  faces  of  the  filters  were 
separated  from  the  image  by  tilting  these  surfaces  slightly.  Lenses,  (Thorlabs  CllOTM-A) 
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with  a  6.24  mm  focal  length  and  a  0.40  numerical  aperture,  focus  each  beam  to  a 
diffraction  limited  spot  with  a  radius  of  about  5  |im  on  the  object  and  reference  mirror. 
These  lenses  are  AR  coated  to  pass  wavelengths  between  350-600  run.  Again,  the  reason 
a  lens  is  used  in  the  reference  beam  is  to  compensate  for  the  dispersion  which  would  be 
caused  by  a  single  lens  in  the  object  beam  The  scattered  light  from  the  object  and 
reference  mirror  is  then  re-coUimated  by  their  respective  lenses  and  sent  back  through  the 
beamsplitter  to  be  recombined,  forming  an  interference  pattern  that  can  be  detected  on  a 
CCD  array.  For  practical  reasons,  the  collection  of  the  scattered  light  from  the  object 
differs  from  that  modeled  in  Ref.  [1]  in  that  here  we  are  probing  the  object  in  a  reflective 
geometry  instead  of  a  transmissive  one. 

The  reference  mirror  position  or  path  delay  of  the  reference  beam  is  controlled  by 
two  systems.  First,  we  use  an  electro-motive  translator  (MPB)  with  analog  closed  loop 
feedback  from  an  inductive  position  sensor  (Kaman  Instrumentation)  which  gives 
maximum  stability  and  linearity.  This  system  is  in  turn  controlled  using  a  novel 
interferometric  positioning  system  designed  specifically  for  this  experiment.  This  system 
uses  a  second  laser  (Spectra  Physics  frequency  stabilized  He:Ne  laser)  which  propagates 
slightly  above  and  parallel  to  the  dye-source  beam,  thus  traversing  the  identical  optical 
path  distance.  The  spatial  fringe  pattern  produced  by  this  He:Ne  laser  is  captured  by  the 
CCD  array  and  Fourier  transformed  to  directly  extract  the  phase  of  the  spatial  fringes.  By 
maintaining  a  careful  history  of  this  phase,  the  absolute  position  of  the  reference  mirror 
can  be  found  to  within  1  nm  over  the  full  range  of  the  electro-motive  translator  (100  |J,m.) 
This  positioning  mformation  is  then  fed  back  to  the  analog  feedback  system  with  a  16-bit 
digital  to  analog  converter  to  con:q)lete  the  digital  feedback  loop.  This  results  in  our 
ability  to  determine  the  OPD  to  an  absolute  accuracy  of  less  than  2  nm  and  to  control  the 
OPD  in  steps  of  5  nm.  We  obtain  this  accuracy  over  the  system’s  con^lete  OPD  range  of 
200  |xm  which,  incidentally,  is  twice  the  range  of  the  actual  stage  movement  and  is  due  to 
the  nature  of  the  Michelson  type  interferometer  used. 

A  scientific-grade  CCD  array  (Princeton  Instruments)  is  used  to  capture  the  data 
from  the  cross-correlations  as  well  as  to  capture  the  fringe  pattern  produced  by  the  He:Ne 
laser  for  the  positioning  system  This  CCD  array  features  active  cooling,  16-bits  of 
dynamic  range,  a  full  well  capacity  of  500,000  electrons,  and  an  integrated  electro¬ 
mechanical  shutter.  To  eliminate  excessive  mechanical  vibrations,  we  externally  mount  the 
shutter  to  an  isolated  structure. 

To  laterally  position  the  sanq)le  in  the  x  and  y  directions,  we  use  a  pair  of 
rnmpniind  flexure  Stages  with  piezoelectric  drives  (Melles  Griot  NanoFlex™ 
17TFC(X)1/MD)  which  feature  zero  friction  and  stiction  over  a  20  |im  range.  These 
stages  also  rely  on  analog  closed  loop  feedback  by  employing  inductive  sensors  (Kaman 
Instrumentation)  for  maximum  stability  and  linearity.  The  sensors  are  factory  calibrated 
using  interferometric  techniques  which  give  us  an  accuracy  better  than  50  nm,  assuming 
the  temperature  is  stable  over  the  length  of  the  measurement.  Note  that  only  x  positioning 
is  needed  in  the  experiments  discussed  this  section. 

The  subwavelength  featured  object  used  is  the  semi-polished  back  of  a  GaAs 
wafer.  When  viewed  under  a  scanning  electron  microscope,  as  shown  in  Fig.  3,  the 
surface  is  observed  to  contain  a  random  pattern  of  grooves  and  pits,  a  result  of  the 
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chemical  etch  used  in  processing.  The  sharp  edges  have  high  spatial  frequencies,  which 
provide  adequate  form  dispersion  for  our  system  to  detect. 

This  system  is  very  sensitive  to  noise  and  changes  in  the  surrounding  environment. 
Therefore  everything  possible  must  be  done  to  isolated  the  system  and  reduce  sources  of 
noise.  The  system  is  mounted  on  a  vibration  danced  optical  table  (TMC)  supported  on 
air  pistons.  Enclosures  are  used  to  isolate  the  dye  source  from  the  interferometer  as  well 
as  isolating  the  CCD  array.  The  enclosure  used  for  the  interferometer  is  acoustically  and 
thermally  insulated.  The  temperature  of  this  enclosure  is  held  constant  at  77.26  ®F  using  a 
wire  heater  with  temperature  feedback  control.  All  remaining  optical  paths  are  covered  to 
prevent  dust  and  air  turbulence  caused  by  the  air  handlers  from  becoming  sources  of  noise. 
All  mirrors  used  in  the  system  with  exception  of  the  dye-source  have  a  surface  flatness  to 
A/20  to  insure  minimal  wavefront  distortion.  A  single  con^uter  is  used  to  control  the 
positioning  of  the  reference  mirror  as  well  as  store  data  collected  from  the  CCD  array 
during  the  cross  correlations.  The  computer  is  networked  to  allow  monitoring  from  a 
remote  location  so  as  not  to  disturb  the  experiment  while  taking  data.  A  schematic  of  the 
overall  system  is  shown  in  Fig.  1. 

Our  experimental  procedure  consists  of  collecting  cross-correlations  over  a  square 
aperture  on  the  CCD  array  measuring  11x11  pixels  to  give  a  total  of  121  spati^ 
channels.  The  cross-correlation  of  an  object  as  performed  with  our  interferometric 
imaging  system  is  accomplished  by  indexing  the  OPD  in  steps  over  the  foil  extent  of  the 
“terrporal”  fringes.  The  step  size,  which  is  150  nm,  is  chosen  so  that  the  sampling  of  the 
interference  fringes  is  near  the  Nyquist  rate.  The  extent  of  “ten^Joral  fringes  in  our 
system  will  always  be  greater  than  one  over  the  bandwidth  of  the  source  and  is  found  by 
experin^nt  to  range  between  300  and  1500  samples.  For  each  value  of  the  reference 
mirror  position  or  OPD,  the  intensity  at  each  pfacel  in  the  1 1  x  1 1  aperture  is  integrated  for 
a  specified  timp.  and  then  recorded  by  the  computer.  This  time  duration,  which  is 
controlled  by  a  mechanical  shutter,  ranges  between  50  and  100  ms,  depending  on  the 
incident  intensity,  and  is  optimized  at  the  beginning  of  an  experiment  to  utilize  the  foil  well 
capacity  of  the  CCD  array.  A  sanqile  of  the  interference  fiinges  collected  at  one  pixel 
versus  the  OPD  is  shown  in  Fig.  2. 

Several  cross-correlations  are  collected  for  each  of  several  object  positions  as 
foUows.  We  begin  by  laterally  moving  the  object  to  position  1  and  collecting  a  cross- 
correlation.  The  object  is  then  laterally  shifted  by  Ax  (where  Ax  is  less  than  one-half  of  the 
center  wavelength  of  the  scanning  light)  to  position  2  and  another  cross-correlation  is 
captured.  Cross-correlations  are  collected  for  up  to  five  different  object  positions.  The 
entire  process  is  repeated  up  to  50  times  so  that  redundant  data  is  taken  for  each  object 
position.  This  means  that  each  individual  object  position  can  have  up  to  50  cross- 
correlations  associated  with  it.  To  maximize  the  SNR,  it  is  necessary  to  coUect  as  much 
information  as  is  reasonably  possible.  Two  factors  which  limit  the  amount  of  data  that  can 
be  collected  are,  the  amount  of  data  storage  space  available,  and  the  time  it  takes  to 
coUect  the  It  can  take  up  to  12  minutes  to  collect  a  single  cross-correlation.  Thus 
the  time  it  takes  to  collect  a  foil  set  of  data  (50  cross-correlations)  can  be  as  long  as  10 
hours.  It  is  crucial  to  make  this  time  as  small  as  possible  to  avoid  any  time  dependent 
correlations  from  entering  into  our  data. 
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in.4.  EXPERIMENTAL  RESULTS 

The  following  experimental  results  indicate  subwavelength  detection  of  an  object  using 
conventional  optics  coupled  to  an  interferometric  detection  system  The  resolution  limit  is 
being  tested  with  subwavelength  translation  of  the  object  and  indicates  resolution  to  better 
than  A/10.  Since  the  accuracy  of  the  positioning  stage  that  moves  the  object  begins  to 
drop  off  around  50  nm,  we  are  quite  sure  this  is  affecting  the  lower  bound  on  the  system’s 
resolution.  The  actual  resolution  of  our  system  may  be  substantially  better  than  our  results 
indicate. 

The  results  of  a  particular  experiment  described  in  Section  3  are  shown  in  Fig.  4. 
In  this  experiment,  we  collect  nine  cross-correlations  at  each  of  five  different  object 
positions  separated  by  50  rnn,  for  a  total  of  45  cross-correlations.  The  position  of  the 
object  is  changed  after  each  cross-correlation  was  collected  to  avoid  correlations  between 
the  object  position  and  various  time  dependent  effects  (i.e.  teiiq)erature).  The  data  is  then 
processed  as  described  in  Section  2.  First  a  set  of  averaged  basis  vectors  (j3)  are  calculated 
for  each  object  position.  Then  one  additional  cross-correlation  is  collected  fi'om  each 
object  position  and  compared  to  the  average  basis  vectors  already  computed.  When 
Eq.  (3)  is  applied,  the  a’s  are  taken  fi-om  this  latest  set  of  five  cross-correlations  and  the 
p’s  are  taken  fi’om  the  first  45  cross-correlations.  We  use  separate  sets  of  cross- 
correlations  to  calculate  the  a’s  and  the  ps.  This  is  to  ensure  the  averages  (ps)  are  not 
biased  by  the  cross-correlations  we  want  to  identify.  What  we  are  essentially  doing  by  this 
process,  is  building  a  map  of  nine  different  transfer  functions  of  the  object  at  each  position. 
Then,  by  taking  additional  cross-correlations  at  previous  object  positions  and  determining 
their  transfer  fonctions,  we  can  identify  fiom  which  object  position  these  latest  cross¬ 
correlations  were  taken.  We  have  found  that  for  50  nm  movements  of  the  object  (less 
than  one-tenth  the  wavelength  of  the  scanning  light),  that  we  can  identify  with  100% 
accuracy  the  position  of  the  object.  Rgure  4  shows  the  results  of  this  experiment  for  e^h 

of  five  different  object  positions.  Plots  4(A)  through  4(E)  represent  {7,  ^  ^  through 

’  respectively.  The  minimum  in  each  plot  corresponds  to  the  average  position 

that  is  closest  to  the  actual  position  measured.  The  values  of  in  these  plots  is  5. 
Fourteen  evenly-spaced  fi»quency  continents  between  490.0  THz  and  555.0  THz  are 
used  for  the  analysis.  The  spacing  between  the  fiequencies  coiiinents  is  5  THz.  This 
corresponds  to  a  wavelength  range  of  540  nm  to  612  nm 

Figure  5  shows  the  accuracy  of  our  interferometric  imaging  system  for  different 
step  sizes.  In  each  case,  ten  cross-correlations  fix)m  five  object  positions  are  measured. 
Then,  the  basis  vector  coefficients  used  to  represent  each  of  the  50  transfer  functions  are 
con:q)ared  to  the  average  coefficients  for  each  object  position  as  in  Eq.  (3).  The  percent 
accuracy  shows  the  percentage  of  times  the  object  position,  determined  by  the  minimum 
value  of  7,  correctly  corresponds  to  the  actual  object  position  fiom  which  the  cross¬ 
correlation  is  taken.  The  preceding  analysis  is  done  using  five  basis  vectors  ( =  5 )  and 
14  adjacent  frequency  con^onents. 
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Figure  6  shows  how  the  accuracy  is  affected  when  the  number  of  basis  vectors 
^max  ^  decreased.  This  experiment,  as  above,  had  five  object  positions,  and  ten 
measurements  at  each  position.  The  step  size  was  50  nm  and  the  number  of  adjacent 
frequency  components  is  14.  As  can  be  seen,  the  results  are  correct  for  four  or  more  basis 
vectors  (n^  =4.) 

Figure  7  shows  that  we  can  get  results  with  100%  accuracy  with  as  few  as  8  spatial 
channels,  or  pixels.  In  this  plot,  =  5,  the  step  size  is  50  nm,  and  the  number  of 
adjacent  frequency  conqionents  is  14.  Note  that  these  8  pixels  are  all  in  a  line, 
perpendicular  to  the  direction  of  object  movement.  However,  when  the  pixels  are 
arranged  horizontal  to  the  direction  of  travel,  the  accuracy  decreases  substantially. 

Figure  8  shows  how  the  accuracy  scales  with  the  number  of  fi’equency  components 
used  in  the  analysis.  Since  the  whole  of  the  analysis  is  dependent  on  multiple  fi'equencies, 
we  obviously  should  not  expect  perfect  results  when  using  only  one  fi’equency  con^onent. 
From  this  plot,  we  can  see  that  below  4  frequency  components,  the  results  tend  to  be  less 
than  100%  accurate,  which  confirms  our  theoretical  work.  The  frequency  con^onents  are 
centered  around  the  maximum  of  510.0  THz  and  are  spaced  5.0  THz  apart. 
in.5.  CONCLUSION 

We  have  shown  that  interferometrically  detected  polychromatic  fields  scattered  by  small 
objects  formed  of  nondispersive  materials  contain  sufficient  degrees  of  freedom  for 
imaging  resolution  superior  to  that  of  monochromatic  systems.  In  addition,  we  have 
shown  that  this  extra  information  can  be  detected  experimentally  and  used  for 
superresolved  object  discrimination.  Additional  studies  extending  our  theoretical  analysis 
to  2-D  structures  formed  of  dispersive  scatterers  and  considering  the  problem  of  image 
reconstruction  must  be  performed  to  apply  these  degrees  of  freedom  to  superresolved 
optical  imaging. 

in.  6.FIGURES 

Fig.  1.  Interferometric  imaging  system 

Fig.  2  Intensity  at  one  spatial  location  versus  the  position  of  the  reference  mirror. 

Fig.  3  Semi-poUshed  GaAs  surface  viewed  under  a  scanning  electron  microscope  at 
1600  X  magnification. 

Fig.  4  Y  versus  position  for  five  different  positions.  The  impulse  responses  used  to 
generate  Plots  (A)  through  (E)  were  taken  at  positions  1  through  5, 
respectively.  Since  the  minimum  in  each  plot  corresponds  to  the  position  the 
impnkft  response  was  taken  at,  that  position  was  correctly  determined.  For 
these  plots,  =5,  the  number  of  frequencies  was  14,  and  the  number  of 

pixels  was  121. 

Fig.  5  Accuracy  versus  step  size,  =  5 ,  number  of  fi^quencies  is  14,  and  number 
of  pixels  is  121. 

Fig.  6  Accuracy  versus  number  of  basis  vectors  .  Step  size  is  50  nm,  number  of 
fi:equencies  is  14,  number  of  pixels  is  121. 

Fig.  7  Accuracy  versus  number  of  pixels,  =  5 ,  step  size  is  50  rnn,  number  of 
frequencies  is  14. 
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Fig.  8 


Accuracy  versus  number  of  frequencies. 
of  pixels  is  121. 


=  5 ,  step  size  is  50  nm,  number 
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Fig.  1.  Interferometric  imaging  system. 
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Fig.  3.  Semi-polished  GaAs  surface  viewed  under  a  scanning  electron  microscope 
1600  X  magnification. 


Fig.  4.  7  versus  position  for  five  different  positions.  The  impulse  responses  used  to  generate  Plots  (A) 
through  (E)  were  taken  at  positions  1  through  5,  respectively.  Since  the  minimum  in  each  plot 
corresponds  to  the  position  the  impulse  response  was  taken  at,  that  position  was  correctly  determined.  For 
these  plots,  n,M«  =  5,  the  number  of  frequencies  was  14,  and  the  number  of  pixels  was  121. 
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step  size  (nm) 

Fig.  5.  Accuracy  versus  step  size,  =  5,  number  of  frequencies  =  14,  number  of  pixels  =  121 . 


Fig.  6.  Accuracy  versus  number  of  basis  vectors  rima.  St^  size  =  50  nm,  number  of  frequencies  =  14, 
number  of  pixels  =  121. 
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number  of  pixels 

Fig.  7.  Accuracy  versus  number  of  pixels,  ttm^  -  5,  step  size  =  50  nm,  number  of  frequencies  =  14. 


number  of  frequenciess 

Fig.  8.  Accuracy  versus  number  of  frequencies,  ttmu  =  5,  step  size  =  50  nm,  number  of  pixels  =  121. 
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